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Abstract 

A hierarchical froth model of the interface of a random (/-state Potts ferro- 
magnet in 2D is studied by recursive methods. A fraction p of the nearest 
neighbour bonds is made inaccessible to domain walls by infinitely strong ferro- 
magnetic couplings. Energetic and geometric scaling properties of the interface 
are controlled by zero temperature fixed distributions. For p < p c , the directed 
percolation threshold, the interface behaves as for p = 0, and scaling supports 
random Ising (q = 2) critical behavior for all q's. At p = p c three regimes are 
obtained for different ratios of ferro vs. antiferromagnetic couplings. With rates 
above a threshold value the interface is linear ( fractal dimension df = 1) and 
its energy fluctuations, AE scale with length as AE oc L w , with u — 0.48. 
When the threshold is reached the interface branches at all scales and is fractal 
(df ~ 1.046) with lu c ~ 0.51. Thus, at p c , dilution modifies both low temper- 
ature interfacial properties and critical scaling. Below threshold the interface 
becomes a probe of the backbone geometry (df — d — 1.305; d = backbone 
fractal dimension ), which even controls energy fluctuations (u — df — d). 
Numerical determinations of directed percolation exponents on diamond hierar- 
chical lattice are also presented. 
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1 Introduction 

The effects of quenched disorder due to random impurities on phase transitions 
are often non trivial and represent since long an active research field. According 
to a simple Harris criterion [1] , disorder is relevant to continuous phase transitions 
with positive specific heat exponent. The presence of impurities can also change 
first order transitions into second order, and typically has rounding effects upon 
them [2,3]. Numerical and experimental evidence exists [4,5] that in some random 
systems, like Potts ferromagnets, the critical behaviour falls into the random bond 
Ising model (RBIM) universality class, irrespective of the different symmetry they 
possess. Recently a justification of this possible universality has been proposed, 
based on the analysis of the behaviour of critical interfaces in presence of strong 
disorder [6]. 

The interfacial free energy is indeed crucial to second order phase transitions since 
it vanishes in a singular way at the critical point. RBIM interfaces have been suc- 
cessfully investigated within the solid-on-solid (SOS) model, which ignores over- 
hangs and islands near the interface [7]. This model is equivalent to the directed 
polymer in random medium (DPRM) [8], which reduces to an optimization prob- 
lem of path energies, at T = 0. 

The (/-state Potts model generalizes the Ising one with its S q symmetry. When 
the model undergoes continuous phase transitions, Potts interfaces are compli- 
cated branched objects which become fractal at criticality, in the pure case. The 
so called "froth model" has been introduced as a simplified description of those 
interfaces [6]. This model represents a nontrivial generalization of the DPRM 
and provides a possible mechanism for the explanation of the puzzling universal- 
ity in the strongly disordered Potts model at criticality. The realization of the 
froth model on a diamond hierarchical lattice (DHL) in presence of coupling con- 
stant randomness is the only context where the hard computational difficulties 
connected with its study could be successfully faced so far. Indeed, on DHL the 
model is amenable to an accurate position space renormalization group analysis 
yielding deep insight into interfacial scaling properties [6]. 

A variant to the DPRM has been recently obtained through the introduction of a 
geometrical disorder represented by bond dilution [9] , which shows to be itself an 
implicit source of randomness for path energy distributions. Dilution, conceived as 
the presence of a fraction of bonds inaccessible by the polymer, modifies or leaves 
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unaffected DPRM critical behaviour depending on whether the concentration of 
present (accessible) bonds is at or above the directed percolation threshold, respec- 
tively. The structure of the underlying directed bond percolation cluster governs 
the scaling behaviour of the model right at the critical threshold [9]. 
Interest in investigating the effects of a similar dilution on the froth interface model 
of ref.[6] is suggested by different considerations. As the DPRM can describe frac- 
ture [10], the froth interface can also be seen as realizing a more complicated 
pattern of fracture, with the possibility of branchings and loops formation in the 
cracks. Dilution can thus represent in both cases the effects of hard inclusions 
through which cracks can not propagate, in an inhomogeneous material. In some 
systems like 3 He- 4 He mixtures in the pores of aerogels, dense coatings of 4 He are 
formed on the surface of the aerogel structure [11]. The higher density considerably 
enhances the interaction between He atoms. This effect adds to a more general 
introduction of disorder in the interactions. One expects that enhanced exchange 
interactions around the aerogel prevent penetration of the interface associated to 
the superfluid transition in the aerogel [12]. The space regions occupied by the 
aerogel and its coating can again be seen as inclusions with very hard coupling for 
the interface. 

These examples suggested us to consider a diluted generalization to the froth inter- 
facial model of ref.[6], by allowing lattice bonds to be absent, i.e. unaccessible to 
the interface, with a finite probability. As in ref.[9], we choose here a percolative 
dilution geometry, which is of course rather different from that of the aerogels' 
case. 

Dilution, especially at percolation threshold, is an obvious candidate to possibly 
introduce new universality classes of interfacial scaling behaviour. This is most 
important for branching interfaces in random environment, because their study is 
still in a very preliminary stage [6]. 

A further motivation for diluting the froth interface model comes from the legit- 
imate hope that some of its scaling properties at percolation threshold could be 
more or less directly related to the geometry of the critical percolation cluster. 
This article is organized as follows. In the next section we introduce the froth 
model of the Potts interface and set up an iterative renormalization group (RG) 
analysis for its study on hierarchical lattice. We discuss in particular the specific 
features implied by the presence of dilution. In the third section we present our 
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results for several regimes of dilution. The fourth section is devoted to a summary 
and to concluding remarks. 

2 Diluted model of the Potts interface 
The interface between two distinct ordered phases of a 2D random bond (/-state 
Potts ferromagnet at very low temperature is a slightly fluctuating line. Rising 
the temperature to approach its critical value, bubbles of any of the other q — 2 
states may form at the interface; neglecting isolated islands and overhangs, as 
in SOS models of Ising surfaces, we are left to consider an aggregate of bubbles, 
the froth, each bounded by two SOS surfaces. This froth represents the original 
interface. In a random bond Potts problem the interface crosses (breaks) bonds 
with random couplings, thus getting a random energy contribution at each step. 
In addition to random couplings, here we suppose there exist, with probability p, 
steps inaccessible to the interface (dilution). 

Even starting from the above simplified picture, the corresponding interfacial 
model is too hard to analyse. By confining the interface to the bonds of the 
DHL shown in Figure 2.1, we obtain a simplified version of the froth model, where 
a position space RG strategy can be worked out. This model was already studied 
in ref. [6] in the undiluted case (p = 0) , with only random exchange effects. 
The interface partition function, Z n , on the diamond hierarchical lattice at its 
n— th construction stage (longitudinal length L = 2 n ) is calculated iteratively. In 
the pure case (no randomness and no dilution), 

Z n+1 = 2Z 2 n + (q- 2)Zi (2.1) 

where Z = e-? 2J is the Boltzmann weight of a broken bond, and (3 — l/Zc^T. The 
entropic factor q — 2 accounts for the number of different coexisting phases which 
can occupy a bubble. An analysis of the RG flow described by eq.(2.1) shows that 
there is a finite unstable fixed point Z c , besides Z = and Z = oo. The quantity 
/ = lim ln(Z n )/L should be < and coincide with the interfacial line tension 

n — >oo 

for T < T c . So, the region Z < Z c corresponds to the low temperature regime, 
because there Z = exp (L/), where / (/ < 0) is the line tension. For Z > Z c the 
line tension description does not apply, since upon iterating eq.(2.1) Z approaches 
infinity as exp(L 2 / & ), where ft, > is a sort of dense froth free energy which 
describes the system in the high-T phase. The RG analysis provides the critical 
Boltzmann weight, Z c , and the interface critical exponent fi(q), characterizing the 
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vanishing of the interfacial free energy, / ~ \T — T C | M , for any q [13]. Indeed, if 
we put t n = Z n — Z c , we get by construction /(to) = 2 _1 /(ti) = 2 _1 /(2 yT to), 
which leads to the expected scaling of / for t — > 0, with ,u = ln(||^-)/ In 2 = y T . By 
keeping the lengths of the system finite and equal to L and L' = L/2, we also obtain 
the finite size scaling (FSS) version of the previous law: f(to, L) = 2 _1 /(2 yT to, L'). 
Of course, the thermal exponent yx has also the meaning of fractal dimension of 
the critical interface because the average number of broken bonds is obtained by 
differentiating \nZ n with respect to to- 
rn presence of random couplings and bond dilution, Z n is a random variable whose 
distribution function, P n (Z), has to be iterated starting from the initial form: 

P (Z)=p6(Z) + (l-p)V (Z) . (2.2) 

Namely each accessible bond is given a weight, Z, according to Vo(Z). p is the 
probability that a given bond is not accessible to the interface. For a given real- 
ization of disorder (2.2) in the whole structure, the partition function evolves from 
a stage to the next one as: 
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Z n+1 = Z n (l)Z n (2)+Z n (3)Z n (4) + (q-2)l[Z n (i) , (2.3) 

i=i 

where the indices from 1 to 4 refer to the four n-order elements constituting an 
n + 1-order diamond. P n (Z) iterates according to the rules of the nonlinear com- 
position of independent random variables defined in (2.3). Fixed points of the 
mapping (2.1) are replaced by fixed distributions. 

It is easy to see that the component of P n concerning the probability for a bond 
to be infinitely hard decouples from the other quantities and iterates separately 
[9] as: 

p n+1 =[l-(l-p n ) 2 ] 2 , (2.4) 

where p n is the probability for inaccessible macrobonds to occur through the 
n— level lattice. So (2.4) is nothing but the RG transformation for percolation 
on DHL, whose critical fixed point is at p c /2 [9]. As we will also 

discuss at the end of the next section, this is a reasonably good qualitative model 
of directed percolation in 2D. 

The second component, Vq(Z), of Pq{Z) is related to energy randomness, and 
couples to bond dilution under RG flow. Hence even in non random energy cases, 
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Vq{Z) = 5{Z — Z ), Z\ no longer possesses a binary distribution with values 
and Z, already after the first iteration. Thus P n (Z) always evolves towards a 
multivalued distribution function, like in the random energy model. 
In our model we verified that both dilution and disorder are relevant perturba- 
tions of the pure system fixed point of eq.(2.1). Furthermore, similarly to what 
happens in the undiluted case [6], the impure system behaviours turn out to be 
always controlled by T = fixed point distribution functions. This is similar to 
what happens in analogous problems [8,14], where randomness is relevant and en- 
ergy optimization criteria prevail with respect to entropic considerations. Dilution, 
here, is a further source of disorder. So it is no surprise that the dilute random 
system is governed by zero temperature strong disorder fixed-point distributions 
for all dilution regimes. We examine properties of these distributions in detail 
below by considering from the start T = RG recursions. 

We are then concerned with the iteration of the interfacial energy distribution 
starting from one which introduces absent (i.e. inaccessible) bonds with proba- 
bility p and present (i.e. accessible) ones with probability 1 — p. The energies 
e of the latter at the n = level are chosen according to a probability density: 
V(e) = u8(e — 1) + (1 — u)5(e + 1). In the pure system we had e = 2 J for all bonds. 
By decreasing u we expect to reach criticality conditions similar to those found by 
rising T at nonzero temperature [6] . On a given sample, the interface at level n + 1 
can cross either of the two lattice branches made by two linked n— macrobonds, 
or form a bubble through both. The choice among these three possibilities is only 
dictated by the criterion of energy minimization, provided the n— macrobonds are 
globally accessible. Monte Carlo iteration of the energy distribution density is 
relatively simple and fast [15]. Starting with a large sample (~ 2 ■ 10 5 ) of bond 
energies we could iterate it up to n ~ 20 without propagating too large statistical 
uncertainties due to finite samplings of the distributions. 

3 Numerical and exact results on RG flows 
RG flows confirm the existence of a critical percolation threshold for absent bonds 
at p = p c . This is indeed the threshold above which the interface can not anymore 
cross the system, unless at the cost of an infinite positive energy, even for finite 
system sizes. The numerical value p c = 0.3820 ± 0.0008 compares quite well with 
the theoretical one. As in ref. [9] we find the influence of dilution on the interfacial 
critical system to be qualitatively different depending whether p < p c , or p = p c . 
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In the p < p c regime forbidden regions do not percolate and only isolated finite 
clusters of inaccessible bonds exist. Clusters are fractal on scales shorter than 
their linear size. On large length scales they are homogeneously distributed and 
should not exert any particular influence on the interface, except as an additional 
source of disorder. We observe in the interface two distinct phases separated by a 
critical value, u c , of the disorder parameter u (u/ (1 — u) gives the ratio of ferro- 
vs. antiferromagnetic accessible bonds). u c of course depends on p and on the fact 
that we were setting to unity the starting absolute e value. 

When u > u c the interface is linear, its energy distribution has mean value, (E(L)), 
which shifts as (E(L)) ~ (E )L > (L = 2 n is the length of the DHL at level 
n) and a width, AE(L), growing as L u , with u 0.30. The interface is thus a 
directed path on a DHL with the energy distribution discussed in [16]. For u < u c 
we find a dense foam of bubbles, whose energy density function has mean value 
(E(L)) oc —L 2 and width AE(L) oc L. The central limit theorem applies to this 
distribution function which should be asymptotically Gaussian. At u = u c we 
have (E(L)} oc L Vc and AE(L) oc L u % with y c ~ u c « 0.32, suggesting that on 
approaching u c the average interface energy should scale as: 

(E(p,u,L)} = A\u-u c \^L + BL^ . (3.1) 

Our numerical estimate is \x = 0.93 ± 0.08. By further considering that the crit- 
ical interface is still linear, since its mass, M(L), scales as M(L) oc L, we find 
here remarkable coincidences with the undiluted random bond Potts interface [6] , 
suggesting that the observed transition should fall in the same universality class. 
Indeed, the linearity of the interface, and the value of [i compatible with unity 
suggests that the randomness produces Ising critical behavior independent of q as 
discussed in ref.[6]. This conclusion is corroborated by a study of the probability 
Pi(n) that the interface forms a loop at stage n. Like in the undiluted case, this 
probability tends to as n approaches infinity. 

In the critical percolating regime, p = p c , both absent and present bonds percolate 
and form incipient macroscopically spanning clusters which have strong effects on 
the interface shape. Again two phases exist, separated by a critical value, w c , of 
u. 

The interface is linear in the region above u c and the energy distribution function 
flows under RG to a phase characterized by a mean (E(L)} oc L and a width 
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AE(L) oc with uj = 0.48 ± 0.02, peculiar of a directed path confined on an 
essentially one-dimensional structure (the percolation critical backbone). This is 
precisely the regime found for the linear DPRM at the directed percolation thresh- 
old in ref.[9]. 

When u < u c a dense aggregate of bubbles forms, saturating all paths which 
connect the ends of the DHL, hence invading the whole percolation backbone of 
present bonds. The energy has a distribution whose mean and width scale as L y 
and L u , respectively, with y = uj = 1.305 ± 0.006. Since the backbone has fractal 

dimension d = 1.30575 [17], we argue that the critical percolation geometry 

strongly controls the interface energy fluctuations. Indeed these are of the same 
order as those of the percolation cluster mass (uj ~ d). The fact that the backbone 
mass and its fluctuations alone determine the interface energy scaling properties 
is confirmed by the results for u = (all bond energies equal and negative), which 
yield the same uj ~ d. Thus, backbone mass fluctuations scale like the average 
mass, and the randomness in bond energies does not add new features to the in- 
terface energy fluctuations. 

At the transition point, u = u c = 0.661 ± 0.001, the scaling exponents of (E(L)) 
and AE(L) versus L are, respectively: y c = 0.50 ± 0.02 and ui c = 0.51 ± 0.01. On 
approaching u c the average interfacial energy fulfills a scaling of the form: 

(E(p,u,L)} = A\u-u c \»L + BL^ . (3.2) 

with [i = 0.94 ±0.04, as illustrated by the collapse of data in Figure 3.1. This value 
is not incompatible with the undiluted froth model one [6]. However, collapse fits 
are not the only way to estimate fx. In order to deeper investigate the nature of this 
transition it is worthwhile calculating also the scalings of interface mass, M, and 
probability of forming loops , Pi(n), at level n [6]. M grows as L df . df indicates 
the interface fractal dimension, and we find in the three different regimes: 

u>u c d f = 1.000006 ±3- 10" 6 , 

u = u c d f = 1.046 ±0.006 , 

u<u c df = 1.3054±6- 10" 4 . 

Notice that, as already mentioned, df ~ uj ~ d in the saturated phase, u < u c (see 
also Figure 3.2). 
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The loop probability is asymptotically independent of 7i, Pi — > Poo, and -Pc» — 
2 df ~ 1 — 1 [6] provides an alternative way to estimate df. The following behaviours 
are found under iteration (L = 2 n ): 

u > w c , Pi(n) ~ A ie " 6L , 

which implies P^ = 0, hence = 1. Notice that A\ and 6 depend of course on 
the actual u value chosen. 

u = u c , Pi(ri)~P 00 + — , 

with 7 = 0.40 ± 0.01, P^ = (2.3 ± 0.6) • 10" 2 , thus d f = 1.03 ± 0.02; 

u < u c , Pi(n) ~ Poo - A 3 e" cL , 

with Poo = 0.236 ± 0.0009, which implies d f = 1.306 ± 0.001. Unlike P^, A 3 and 
c show a dependence on it. 

Obviously there is a linear interface for u > u c , while the interface occupies the 
whole percolation backbone when u < u c . Indeed df ~ d. Most interesting is the 
fact that the interface is no longer linear at the critical point u = u c , because df is 
significantly different from 1, in that case, and the loop formation probability does 
not behave as a marginal field, as it happened in the undiluted case [6]. Since p 
is a relevant scaling field, logarithmic corrections existing in the undiluted model 
disappear, and the critical interface is a branched object, thus indicating that we 
are in presence of a different, new type of criticality. 

Moreover the RG mappings of P n {Z) at finite temperature evolves towards T = 
distributions both within stable phase regions of parameter space, and on ap- 
proaching the critical surface, showing that these phase transitions are also gov- 
erned by zero temperature strong disorder fixed points, as anticipated. 
For p > p c there are of course no paths joining the ends of the DHL, because 
accessible bonds do not percolate anymore. 

Let us now discuss in more detail the percolative landscape across which the inter- 
face develops when p = p c . The probability of a present bond is g = 1 — p, which 
transforms under RG as: 

Qn+l = £>n( 2 - Qn) 2 ■ (3-3) 
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The fixed point g c = 1 — p c sets the percolation threshold. So at p = p c a critical 
cluster of present bonds also exists, as already noticed. Considering the anisotropy 
of DHL, which naturally introduces a preferred direction, the longitudinal one, for 
the percolation process (RG recursions (2.4) and (3.3) are explicit consequences of 
such a directedness), and taking into account the DHL fractal dimension, equal to 
2, the critical percolating cluster on it is eligible to represent a sufficiently faithful 
approximation of the actual structure of the 2D directed percolation cluster. The 
better such a circumstance can be confirmed, of course, the more value should 
be added to our investigation of interfacial properties on DHL. Of course, there 
are limits to such an attitude: it should also not be forgotten, for example, that 
isotropic and directed percolation are not distinguishable on the DHL. 
Since the Potts interface in its saturated phase occupies the whole directed perco- 
lation backbone, we are able to use such interface as a probe of backbone structure 
by the RG analysis. The directed backbone fractal dimension, d, has already been 
determined both analytically [17] and numerically by evaluating the saturated in- 
terface mass scaling, shown in Figure 3.2. We found a value surprisingly close to 
the one known for two-dimensional directed percolation (d ~ 1.31 [18]). 
The correlation exponents of percolation, u n and v ± , are equally accessible on 
DHL. From (3.3) we calculate the rescaling of Ag = g — g c : 

Ag' = X c Ag = b Ve Ag ; A c = 6 - 2Vb . (3.4) 

The two correlation lengths, parallel and normal to the preferred direction, diverge 
near g c as: 

^~\q-q c \-^ ; U~\q-Qc\- v± . (3.5) 

Manifestly, on this hierarchical lattice £y concerns the correlation in the longitu- 
dinal direction connecting the two ends of the lattice. Thus, if £jj = 6 _1 £n, with 
6 = 2 (length rescaling factor), it follows: 

— = = 0.61151617... and v» = 1.6352797... , 

z^i In b 

(the two-dimensional exponent is slightly larger, z/y = 1.735 [19]). 

To derive u ± we define the distance along the direction orthogonal to the lattice, 

as equal to the number of branching levels separating two given sites. With such a 
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convention the average transverse size of the percolation backbone, L ± , transforms 
then according to: 

L(r+1) = 2^(1 - g n ) 2 L^ + 2qI(1 - g n )3L™ + g 4 n 2L^ _ 

0n+l 

Applying finite size scaling to (3.5) we get £ x ~ L v ^l v w, at g = g c . From (3.6) 

, then — = — — — = — , and v ± = 1 . These 
I'll In 2 z/|| 

z/ x 

results should be compared with the two-dimensional exponents: — = 0.633; 

^ii 

u ± = 1.097 [19]. 

Finally the simplicity of our RG also allows for iterating the number of steps on 
which the saturated interface is linear, and which therefore connect different blobs 
of the directed backbone. From a percolation point of view those are just the "red 
bonds" forming one dimensional links of the directed backbone. Like in isotropic 
percolation [20] one can show that the directed backbone red bond mass, Mr, at 

g c has a scaling form Mr oc L Dr , with D R = — , when the system has linear size 

^ii 

L. We measured a numerical value D R = 0.6111 ± 0.0006, which fairly agrees with 
the above theoretic prediction of — . Figure 3.3 shows a plot of Mr versus L. 
The percolation critical exponents calculated above for our simplified model of di- 
rected percolation are remarkably consistent and close to the truly two-dimensional 
ones, that we hope the same could apply to the scaling properties of the Potts in- 
terface we have determined here. 
4 Conclusions 

In this work we addressed the critical behavior of a Potts ferromagnet in the 
presence of both exchange disorder and dilution. The scaling properties of the 
interfacial free energy were considered and dilution acted in the form of percola- 
tive regions not accessible to the interface, due to extremely strong ferromagnetic 
couplings. By restricting interfacial configurations to those which can be hosted 
by a DHL, a very accurate RG analysis of the line tension can be carried on. 
For all dilution regimes considered interface scaling turns out to be controlled by 
T = fixed point distributions. This makes the analysis relatively more easy, since 
T = recursions for the interfacial energy can be studied more effectively than 
their T > counterparts. 

When accessible regions are above their percolation threshold (p < p c ), dilution 
does not seem to modify interfacial scaling with respect to the undiluted case. 
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The three regimes occurring in this case are those on the basis of which the Ising- 
like nature of criticality in disordered Potts ferromagnet in 2D could be argued 
in ref.[6]. New scalings occur when dilution is at threshold and the interface has 
to develop within the incipient infinite cluster backbone. When, at T = 0, u is 
such to guarantee a sufficient dominance of positive energy bonds, the interface 
remains linear at large scales and appears to behave in the same way a strictly 
linear DPRM would in such an environment [9] . Like in the corresponding regime 
above percolation threshold, ramification is fully inhibited if u > u c . However, 
the exponent uj ~ 0.48 falls in a new universality determined by the backbone 
geometry limiting the linear interface [9]. 

A completely new behavior is realized right at u = u c when dilution is at threshold. 
In such conditions we get clear evidence that the interface has nonzero probability 
of branching and thus behaves as a fractal. Its dimension slightly, but definitely 
exceeds 1. Of course, by construction, on DHL a strictly linear interface can not 
behave as a fractal. However, a fractal dimension equal to 1 would persist in the 
linear regime when replacing the DHL with more complicated lattices allowing in 
principle for a fractal linear DPRM geometry. 

The above interface scaling regime at u = u c is the first critical one with nontrivial 
fractal geometry met so far in this kind of studies, and indicates the highly non- 
trivial effect that backbone geometry can exert on the interface and on the critical 
behavior of the system. So far, borderline regimes between the linear and max- 
imally branched ones in such froth models were never seen to allow for nonzero, 
intermediate looping probability. Thus the critical behavior of the interface could 
not be seen to be different from the linear, Ising-like one at low-T [6] . The values 
of \i and df estimated for this regime lead to v = 1/ ji « 1.06 or v = df w 1.046 
for a critical Potts model subject to threshold dilution and disorder. This can be 
concluded within the same limits of the arguments developed in ref.[6]. Like in 
that undiluted case, the value of z/, being associated to a T = fixed distribution, 
is universal with respect to q. 

In the last regime, with critical dilution and u < u c , we could realize how the 
interface can become a sort of probe for the structural properties of the backbone. 
The interface mass M scales as the backbone mass with longitudinal distance. 
Quite remarkably, the energy fluctuations are fully determined by backbone mass 
fluctuations here. These backbone mass fluctuations in turn scale with the same 
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exponent as the total average mass. 

As an instructive exercise adding credibility to our DHL model as a good qualita- 
tive picture of the situation on 2D Euclidean lattice, we produced here results for 
I'll and u ± and the "red bond" dimension of directed percolation. No similar deter- 
mination of u ± was produced before, to our knowledge. These results altogether 
appear remarkably consistent with the numerical estimates in 2D [19]. 
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FIGURES 



FIGURE 2.1. a) The construction rule of the diamond hierarchical lattice. At every stage 
each bond is replaced by a 4-bond cell, b) The DHL at its second stage. 



FIGURE 3.1. Data collapse of (E(u,L,p = p c ))/L Wc against \u - w c | M L 1_Wc . We measure 
H = 0.94 ± 0.04 and uj c = 0.51 ± 0.01. Scales are logarithmic on both axes. 



FIGURE 3.2. Plot of the interface mass, M, versus lattice longitudinal length, L, in the sat- 
urated phase, u < u c , and at critical dilution p = p c . M equals the directed 
percolation backbone mass, which thus exhibits the same scaling against L. 



FIGURE 3.3. Plot of the "red bond" mass, Mr, in the critical percolation backbone, versus 
the lattice longitudinal length L. 
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